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Abstract 



The fourth Standard Model (SM) family quarks and weak iso-singlet quarks 
predicted by Eg GUT are considered. The spin-average of the pseudoscalar 
r]4{n^So) and vector ^4(71^ 5i) quarkonium binding masses of the new mesons 
formed by the fourth Standard Model (SM) family and iso-singlet Eg with 
their mixings to ordinary quarks are investigated. Further, the fine and hy- 
perfine mass splittings of the these states are also calculated. We solved the 
Schrodinger equation with logarithmic and Martin potentials using the Shifted 
large-N expansion technique. Our results are compared with other models to 
gauge the reliability of the predictions and point out differences. 



I. INTRODUCTION 

The toponium quark does not form hadronic states because of its large mass value (m^ ~ 175 
GeV) and full strength of tbW vertex. On the other hand, there are strong reasons that the 
fourth Standard Model (SM) family should exist [1,2]. The flavor democracy (Democratic 
Mass Matrix approach) [2] favors the existence of the nearly degenerate fourth SM family, 
whereas the fifth SM family is disfavored both by the mass phenomenology and precision tests 
of the SM [3] . The fourth SM family fermions and also the isosinglet quarks production at 
colliders have been investigated [3] . Thus, SM may be treated as an effective theory of 
fundamental interactions rather than fundamental particles. The multi-hundreds GeV fourth 
generation up- type quark {u^), if exist, will be produced at the CERN Large Hadron Collider 
(LHC) [3,4] via gluon-gluon fusion [4]. Hence, the observation of the fourth SM family quark 
in ATLAS has been considered in [4,5]. It is expected that the masses of the fourth family 
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quarks lies between 300 < m^^ < 700 GeV with preferable value m^^ = Agwrj — 8m^ = 640 
GeV [3]. The partial- wave unitarity leads to mq < 700GeV = Arut and in general rut <S 
iTT-u^ <^ 777.5. The fourth generation up- type quark would predominatly decay via 774 — > Wb, 
with an expected event topologies are similar to those for quark pair production. The 
best channel for observing it will be [5]: 1*4774 — > WWbb {lu){jj){bb), where one W decays 
leptonically and the other hadronically. The mass resolution is estimated to be 20(40) GeV 
for 777u4 = 320(640) GeV. The pseudoscalar quarkonium state 774(^-S'o) formed by the SM 
fourth generation quarks is the best candidate among the fourth generation quarkonia to be 
produced at the LHC and VLHC [6]. On the other hand, new heavy quarks known as weak 
iso-singlet quarks are also predicted by various extensions of SM and by Eq GUT [7] which 
is favored by superstring theory [8]. For down-type quark (^4), the decay mode is: ^4 — > Wt 
and the final state contains four W bosons: (^4^4 tW~tW+ bW'^W^bW'W'^. 

The small inter-family mixings [9] leads to the formation of the fourth family quarkonia. 
The most promising candidate for LHC is the pseudoscalar quarkonium state, 774, which 
will be produced reasonantly via gluon-gluon fusion particularly through decay channel: 
774 ^ ZH [10]. 

In spite that the masses of new quarks are larger than 777^, they can form new hadrons be- 
cause of the small inter-family mixings leads to the formation of the fourth family quarkonia 
between new heavy and ordinary quarks (t*, d, s, c, b). Indeed, according to the parametriza- 
tion of mass matrices given in [9] mixing between the fourth and third family quarks is 
predicted to be \ Vqu^\ ~ 10^^. Similar situation is expected for iso-singlet down-type quarks. 
The condition for forming {QQ) quarkonia states with new heavy quarks is [11] 

ruQ < {125 GeV) \VQg\-^/\ (1) 

where q — d, s,b for Q — U4 and q — u,c,t for Q — d^. Differing from t quark, fourth family 
quarks will form quarkonia because 7*4 and d^ are almost degenerate and their decays are 
suppressed by small mixings [2,9,12]. 

One of the important goals of the present work is to extend the shifted large-N expansion 
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technique (SLNET) developed for the Schrodinger wave equation [13, 14, 15] and then apphed 
to semirelativistic wave equations [16,17] to reproduce the spectroscopy of the fourth SM 
generation up-type (^4) quark and weak iso-singlet down-type (^4) quark with their small 
mixings to the ordinary type {u, d, s, c, b) quarks. 

Here, we study the present status of the new heavy mesons formed by new quarks in 
the framework of nonrelativistic potential model and give some predictions for their bound 
energy masses. In section II, we present the solution of the Schrodinger equation using the 
flavor-independent logarithmic and Martin potentials for the the self- and non-self conjugate 
Qiq^ and qiqj mass spectra, respectively. A brief conclusion appear in section III. 

II. NEW HEAVY MESONS 

The fourth SM family and Eg isosinglet quarks have formed hadron states if their mix- 
ing with ordinary known quarks is sufficiently small. In the fourth family quarks, the 
parametrization given in [9] satisfies condition (1), whereas new hadrons are not formed 
in the case of paramerization given in [1]. Concerning Eg iso-singlet quarks, we have no 
similar parametrization (one deals with 6x6 mass matrix) and one can make qualitative 
estimations. For example, if the fightest isosinglet quark has the mass rrid^ = 0.5 TeV, new 
heavy hadrons are formed for 1^9^4 1 < 0.09. In order to calculate the binding masses of new 
mesons we investigate two potentials: (1) Logarithmic potential of the form [14,16,18] 

V{r) = -0.6631 + 0.733 ln(r x iGeV), (2) 

where r is the interquark distance which is singular at r = 0. (2) Martin's potential of the 
form [14,16,19] 

V{r) = -8.093 + 6.898r°\ (3) 

which behaves in some respects, as its power approaches to zero, as the logarithmic potential. 
Hence, these potential types are stimulated by the approximate equality of mass splittings 
[18] 
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M{T{bb)) - M{T{bb)) ^ M{tlj'{cc)) - M{tlj{cc)), (4) 

which is independent upon reduced mass. These static quarkonium potentials are monotone 
nondecreasing, and concave functions satisfying the condition [14] 

V'(r) > and V"(r) < 0. (5) 

For two particles system, we shall consider the N— dimensional space Schrodinger equa- 
tion for any spherically symmetric potential V{r). If ip{r) denotes the Schrodinger 's wave 
function, a separation of variables ■^(r) = Y£^m{9,(j))u{r)/r^^~^^^'^ gives the following radial 
equation = c = 1) [13,14,15,16,17] 

r 1 [k - h - a)][k - (3 - a)] ^ 

where /i — {mq^mq.) /{niq^ + ruq.) is the reduced mass for the two quarkonium composite 
particles. Here, En^i denotes the Schrodinger binding energy of meson, and k — N + 2i — a, 
with a representing a proper shift to be calculated later on and / is the angular quantum 
number. We follow the shifted 1/N or 1/k expansion method [13,14,15] by defining 



m=0 

and also the energy eigenvalue expansion [14,15] 

oo -j^(2~m) 

En,l ——Em, (8) 

m=0 ^ 

—1/2 

where x — k (r/ro — 1) with tq is an arbitrary point where the Taylor's expansions is 
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being performed about and Q is a scale parameter to be set equal to k at the end of our 

calculations. Following the approach presented by Ref.[14], we give the necessary expressions 
for calculating the binding energies: 



Q_ 
16/iro 



Eo = V{ro) + (9) 



F - ^ 



1\ (2 -a) 

Ur + - \ 00 — 



2 J 811 



(10) 



Q 



[l-a){3-a) , 



16/1 



(11) 



(12) 



where a*^^^ and a^^^ are two useful expressions given by Imbo et al [13] and also the scale 
parameter Q is defined by the relation 



g = 8/xroV(ro). 



(13) 



Thus, for the A'" = 3 physical space, the Schrodinger binding energy to the third order is 
[14] 



En,e = V{ro) + l:roV'{ro) + \ 



16/x k Vfc 



(14) 



where the shifting parameter, a, is defined by 



a = 2 - {2nr + 1) 
and the root, Tq, is being determined via 



3 + 



roV"{ro) 
V'iro) 



.1/2 



(15) 



1 + 21 + {2nr + 1) 



3 + 



roV"{ro) 
V'iro) 



= [8/.roV'(ro)]'/\ 



(16) 



where = n — 1 is the radial quantum number and n is the principal quantum number. 
Once To is found via equation (16), then the Schrodinger binding energy of the qiqj system 
in (14) becomes relatively simple and straightforward. Hence, the bound state mass of the 
Qigj system is written as 



M{qiqj)ni = mg. + rUg. + 2En,i. 



(17) 



The expansion parameter 1/N or 1/k becomes smaller as I becomes larger since the param- 
eter k is proportional to n which it appears in the denominator in higher-order correction. 

Since systems that we investigate in the present work are often considered as nonrela- 
tivistic system, then our treatment is based upon Schrodinger equation with a Hamiltonian 
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Ho^-^ + V{r) + VsD. (18) 

where Vsd is the spin-dependent term taking the simple form (cf. Ref.[14] and the references 
therein) 

Vsd Vss = ^^^5'{v)s,.s,. (19) 

The spin dependent potential is simply a spin-spin part and this would enable us to make 
some preliminary calculations of the energies of the lowest two S-states of the new quarks 
and their small mixings with the ordinary quarks. The potential parameters in this section 
are all strictly flavor-independents and are fitted to the low-lying energy levels of cc and 66 
systems. The strong coupling constant ag is fitted to the observed charmonium hyperfine 
splitting of 117 MeV. The numerical values of ag for the two potential types have been 
adjusted in Ref.[14] as 

Q;f)(cc) ~ 0.220 and a(^)(cc) ~ 0.251, (20) 

which are found to be dependent on the potential type. Baldicchi and Prosperi [20] used 
the standard running QCD coupling expression 

47r 

(ll-|n,)ln(¥) 

with Uf — A and A = 0.2 GeV cut at a maximum value q;s(0) = 0.35, to give the right 
J/V' — r]c splitting of cc quarkonium and to treat properly the infrared region. Detail on 
their numerical works are given in Ref. [20]. Whereas Brambilla and Vairo [21] took in their 
perturbative analysis 0.26 < as{jji — 2GeV) < 0.30. 

The potential parameters in (2) and (3) together with the quark masses are obtained 
from the experimental SAD mass: M{i/j{lS)) (cf. e.g. Refs.[14,16]). In our calculations 
we used m„ = = 0.367 GeV, rus — 0.561 GeV whereas rric and m;, masses are given by 
Ref. [14]. Further, the fourth SM family up-typc quark and isosinglet down-type quark 
masses are taken as m„4 = 638.6 GeV [9,22] and rrid^ — 0.5 TeV [22], respectively. 



Firstly, we calculate the masses of bound states of cc and bb in Table 1 for the two 
potential types. It is seen that our results are in good agreement with the experimental 
results [23,24]. In Table 2, we present the masses of bound states of quarkonia and mesons 
formed by the fourth SM family up-type quark and Eq isosinglet quark, respectively. In 
Table 3, the predicted mass sphttings is compared with Ref.[22]. In Table 4, the hyperfine 
mass splittings ip4 — r]4 for the first two 5'— states are also given. The best mechanism for the 
production of heavy mesons formed by and quarks is the resonance formation of 35' and 
AS quarkonia at lepton colliders with subsequent decay into corresponding meson-antimeson 
states. 

III. CONCLUSION 

We have produced the mesons formed by new heavy quarks. The mass splitting of the 
logarithmic potentials is found to be independent on the reduced mass and is constant in 
magnitude for any chosen state for all studied quarkonium families. However, it is seen in 
Martin's potential that the mass splittings are equal for any chosen state in the fourth family, 
^4, quark and isosinglet quark with the same type of ordinary quark mixing. Further, the 
hyperfine splitting mass is found to be much smaller compared to the large mass of the new 
mesons. Therefore, the pseudoscalar r]4^{n^So) and the vector ip4{n'^Si) have nearly same 
masses. On the other hand, the Cornell potential fails to produce the masses of the u^u^ 
and (^4^4 quarkonium famileis and their hyperfine mass splittings properly because of the 
small values of the roots ro for the lowest states where the Coulombic part of the potential 
goes to high order fas we take higher order derivatives in the SLNET method. 
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TABLES 

TABLE I. The pseudoscalar, vector and the spin-averaged masses together with hyperfine 

splittings for tlic first two-states of cc and bb states (in MeV). 



States 


cc 


bb 


cc 


bb 


cc 


bb 




Logarithmic: 




Martin: 




Experiment [23,24]: 




IS 


3068 


9444 


3068 


9445 


3068 ± 2 


9448 ± 5 




3097 


9460 


3097 


9461 


3097 


9460 




2980 


9395 


2980 


9397 


2980 


— 




117 


65 


117 


64 


117 


10017 ±5 


2S 


3654 


10030 


3670 


10018 


3663 ± 5 


10023 


2^Si 


3668 


10037 


3685 


10030 


3686 




2'So 


3614 


10008 


3625 


9994 


3622 ± 12^ 




A25 


54 


30 


60 


33 


57 ±8^ 




3S 


3976 


10352 


4021 


10351 




10350 ± 5 


AS 


4200 


10575 


4272 


10590 




10580 


5S 


4370 


10746 


4469 


10777 






IP 


3505 


9881 


3505 


9861 


3525 ± 1 


9900 ± 1 


2P 


3878 


10254 


3907 


10243 




10261 ± 1 



^Here we cite Ref. [25] . 
•^Here we cite Ref. [26]. 



11 



TABLE IL Predicted Spin-averaged masses of the bound states formed by the fourth SM family 
U4 and isosinglet ^4 quarks (in GeV). 



DAL) iviass 


U4U4 




U4U 




1145 




114 c 




1446 




rl .rl . 
04(24 


d4U 




(I4S 




d4C 




rl .h 
CI4O 




Logarithmic: 








































IVl [Ld ) 


LZIO. 




Doy. 


,oU 


Doy. 


,o4 


R'iO 

boy. 


00 

.yz 


RAO 
D4z. 


,oy 


yy / .y4 


oUU, 


. (U 


OUU, 


7/1 
. / 4 


KM 
OUi, 


QO 
.OZ 


Kr\A 
0U4, 


00 

.zy 


M [ZD ) 


LZIO. 


,04 


boy. 


QQ 
,00 


boy. 


,yz 


b4U, 


.OU 


RAQ 

b4o. 


A 1 

,4/ 


one c^Q 
yyo.oo 


OUi, 


OQ 
.ZO 


OUi, 


QO 
.OZ 


OUi, 


on 

.yu 


KCiA 

0U4, 


QQ 
.00 


M{SS) 


1275. 


,96 


640. 


,21 


640. 


,24 


640, 


.82 


643. 


,80 


998.85 


501, 


.61 


501, 


.64 


502, 


.22 


505, 


.20 


M(45) 


1276. 


,18 


640. 


,43 


640. 


,47 


641, 


.05 


644. 


,02 


999.07 


501, 


.83 


501, 


.87 


502, 


.45 


505, 


.42 


iVI l^i-t ) 


LZlO. 


AO 

,4y 


boy. 


, 1 


boy. 


11 
. 1 1 


RAfi 

b4U, 


.00 


RAQ 

b4o. 


QQ 
,00 


nnc QQ 
yyo.oo 


OUi, 


1 Q 

.io 


OUi, 


1 7 
.i 1 


OUi, 


IK 
. 10 


KCiA 

0U4, 


7Q 
. / 


IVI ) 




,OD 


b4U. 


1 1 

,il 


RAr\ 
b4U. 


.10 


b4U, 


. / 


RA Q 

b4o. 


, ( U 


OQQ '7K 

yyo. (0 


OUi, 


.Oi 


OUi, 


KK 
.00 


OUz, 


1 Q 
.io 


OUO, 


1 n 
.iU 


IVI \oI^ ) 


1 07R 


1 1 

,li 


b4U. 


,00 


RAH 
b4U. 


QQ 

,oy 


RAr\ 
b4U, 


.y ( 


RA Q 

b4o. 


,yo 


yyy.UU 


OUi, 


IR 

. ( b 


OUi, 


70 

. / y 


OUz, 


Q7 
.0 / 


OUO, 


IK 
.60 


Martin: 








































IVI[Ld ) 


iz 1 4. 


,oz 


boo. 


, / ( 


boo. 


,oU 


boy, 


.bz 


RAO 

b4z. 


RA 

,b4 


nOT RC\ 

yy /.by 


OUU, 


1 7 
.i t 


OUU, 


on 
.zU 


KM 

oUi, 


no 
.Uz 


Kr\A 
0U4, 


CiA 

.U4 


IVi \Zo ) 




OQ 
,ZO 


boy. 


,oy 


boy. 


,4z 


RAn 
b4U, 


.zi 


RAQ 

b4o. 


on 
,zU 


OOQ 1 ^ 

yyo.io 


OUU, 




OUU, 


QO 
.OZ 


KM 

OUi, 


R'i 

.bi 


0U4, 


Rr\ 
.bU 


M{3S) 


1275. 


,54 


639. 


,76 


639. 


,78 


640, 


.55 


643. 


,52 


998.42 


501, 


.16 


501, 


.18 


501, 


.95 


504, 


.92 


M{4S) 


1275. 


,73 


640. 


,02 


640. 


,03 


640, 


.79 


643. 


,75 


998.61 


501, 


.42 


501, 


.43 


502, 


.19 


505, 


.15 


M(1P) 


1275. 


,15 


639. 


,22 


639. 


,25 


640, 


.05 


643. 


,04 


998.03 


500, 


.62 


500, 


.65 


501, 


.45 


504, 


.44 


M(2P) 


1275. 


,46 


639. 


,64 


639. 


,66 


640, 


.44 


643. 


,41 


998.33 


501, 


.04 


501, 


.06 


501, 


.84 


504, 


.81 


M(3P) 


1275. 


,67 


639. 


,93 


639. 


,94 


640, 


.70 


643. 


,67 


998.55 


501, 


.33 


501, 


.34 


502, 


.11 


505, 


.07 



^Fourth family system scale offset by 2.05 GeV to agree with Ref.[22]. 



'^Iso-singlet system scale offset by 1.96 GeV to agree with Ref.[22]. 
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TABLE III. Predicted fine-splittings of S and P levels (in MeV). 

Mass splittings [22] u^u [22] u^s [22] UiC uj) didi [22] diU [22] ^45 [22] diC dj) 

Logarithmic: 



2S- 


IS 


59 


59 


58 


59 


58 


59 


58 


58 


59 


59 


58 


59 


58 


59 


58 


59 


ZS- 


2S 


32 


32 


33 


33 


32 


32 


32 


33 


32 


32 


33 


33 


32 


32 


32 


32 


AS- 


3S 


22 


22 


22 


22 


23 


23 


23 


22 


22 


22 


22 


22 


23 


23 


23 


22 


2P- 


IP 


37 




38 




38 




38 


37 


37 




38 




38 




38 


37 


3P- 


2P 


25 




25 




24 




24 


25 


25 




25 




24 




24 


25 



Martin: 



2S-1S 


46 


62 


62 


59 


56 


46 


62 


62 


59 


56 


3S-2S 


26 


37 


36 


34 


32 


27 


37 


36 


34 


32 


4S-3S 


19 


26 


25 


24 


23 


20 


26 


25 


24 


23 


2P - IP 


31 


42 


41 


39 


37 


30 


42 


41 


39 


37 


3P-2P 


21 


29 


28 


26 


26 


22 


29 


28 


27 


26 



TABLE IV. Hyperfine mass splittings (in MeV). 



Level 


M4M4 


U4U 


U4S 


M4C 


M46 


did 4 


diU 


d4S 


d4C 


0^46 cc 


hh 


Logarithmic: 


























5.69 


0.39 


0.48 


0.78 


1.39 


6.43 


0.49 


0.61 


0.99 


1.78 117(117)^ 


65(86) 




2.62 


0.18 


0.22 


0.36 


0.64 


2.96 


0.23 


0.28 


0.46 


0.82 54(57 ±8)*^ 


30(35) 


Martin: 


























4.09 


0.45 


0.54 


0.89 


1.38 


4.70 


0.57 


0.69 


1.13 


1.76 117 


64 


A2S 
^The quantity 


2.09 0.23 0.28 
in bracket is the < 


0.45 0.71 2.41 0.29 0.35 0.58 0.90 60 
experimental or taken from other references. 


33 



•^Here we cite Ref.[26]. 
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